Ionization of Rydberg atoms embedded in an ultracold plasma 
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We have studied the behavior of cold Rydberg atoms embedded in an ultracold plasma. We 
demonstrate that even deeply bound Rydberg atoms are completely ionized in such an environment, 
due to electron collisions. Using a fast pulse extraction of the electrons from the plasma we found 
that the number of excess positive charges, which is directly related to the electron temperature T e , 
is not strongly affected by the ionization of the Rydberg atoms. Assuming a Michie-King equilibrium 
distribution, in analogy with globular star cluster dynamics, we estimate T e . Without concluding 
on heating or cooling of the plasma by the Rydberg atoms, we discuss the range for changing the 
plasma temperature by adding Rydberg atoms. 

PACS numbers: 32.80.Pj, 52.25.Dg, 98.10.+Z 



One challenge in ultra-cold plasma physics is to reach 
the correlated regime where the Coulomb energy dom- 
inates the kinetic energy. One suggested way, unfortu- 
nately limited to non alkali ions, is to cool the plasma 
ions by lasers pj. An alternative way might be to use 
the binding energy of Rydberg states as "ice cubes" to 
cool the plasma. The physics of ultracold Rydberg gases 
@>0] and ultracold plasmas Q formed by laser excitation 
of a cold atomic sample have strong similarities. Indeed, 
Rydberg atom formation in an ultracold plasma [j| and 
spontaneous evolution of an ultracold Rydberg gas to 
plasma have been demonstrated. The Rydberg ioniza- 
tion process starts with blackbody photoionization and 
initial electrons leave the cloud region. A second phase 
occurs when the positive ion potential is deep enough to 
trap subsequent electrons, which then collide with Ryd- 
berg atoms creatingmore electrons in an avalanche ion- 
ization process 0,H,EJ- However, other relevant pro- 
cesses have been proposed whose effects need to be in- 
vestigated, such as continuum lowerin g |lfj. and many- 
body effects or long-range interactions IllJ that can lead 
to autoionization of Rydberg atom pairs 12] . 

In this letter, we report the behavior of a mixture of 
an almost neutral ultracold plasma and a cold Rydberg 
atom sample. A related experiment has been reported, 
but in the case of rubidium Rydberg atoms created in a 
purely ionic plasma [l3[ . In this letter, we analyze the 
fast avalanche ionization of deeply bound Rydberg states 
embedded in a quasi-neutral plasma. We also study, with 
a theory based on analogy with globular star cluster dy- 
namics, the evolution of the temperature of the plasma 
when Rydberg atoms are added. 

The cesium magneto-optical trap (MOT) apparatus 
has been described in a previous paper |6| . Two dye lasers 
pulses (Coumarin 500) that are focused to the cold atom 
cloud diameter excite atoms initially in the 6^3/2 state. 
The time origin of the experiment is the first laser (Li) 
pulse, with typical energy P\ — 10 /iJ, which creates a 
quasi-neutral plasma of JVj ~ 4 x 10 5 ions with peak den- 
sity 10 10 cm- 3 . The second laser (L 2 ) pulse (ASE < 1%), 
has a 18 ns delay and excites typically 4 x 10 5 Rydberg 
atoms. The Rydberg number fluctuates from pulse to 



pulse due to the changing overlap of the dye laser mode 
structure with the narrow 6^3/2 — > Rydberg resonance 
|6(. However, the plasma created by the first laser, which 
is tuned just above the ioniziation limit, is affected only 
by negligible laser intensity fluctuations. The MOT trap- 
ping lasers are turned off and a resonant 852 nm diode 
laser pulse excites the Cs 6p 3 / 2 state just before the dye 
laser pulses arrive. Because of the Doppler effect, this 
pulse ensures that all the 6p 3 / 2 atoms are cold 
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have verified that 6s atoms have no effect on our exper- 
imental results by pushing them away using this laser 
light pressure after the pulse dye laser excitation. 

The first part of our analysis concerns the ionizing 
effect of the plasma on the Rydberg atoms. In order 
to study the evolution of the cloud we applied, after a 
variable time t\ =0—20 /its, a positive top hat volt- 
age pulse (voltage V\) to one of the two parallel grids 
(d = 1.57 mm separation) that surround the cloud. This 
pulse is used to pull all the electrons from the plasma 
(hereafter termed "free" electrons) toward a microchan- 
nel plate (MCP) detector, leaving no further electrons 
trapped in the ion space charge in the situation where 
the laser L2 is blocked. The MCP signal is monitored 
using a gated integrator (Gil) (see inset in FigJIJ. The 
effect of the Rydberg atoms on the plasma is monitored 
using a second voltage pulse applied on the opposite grid, 
at time t2 typically 700 ns delay from the first voltage 
pulse. The early part of this second pulse extracts addi- 
tional electrons that are monitored using a second inte- 
grator (GI2). If the second voltage pulse is large enough 
it will also field ionize Rydberg states, leading to a signal 
on the MCP that is monitored using a third integrator, 
GI3. Fig^ shows the results, for t\ — 4 /is, versus the 
wavelength of L2. Lj is blocked by a mechanical shut- 
ter half of the time so that we obtained spectra with, 
and without, the plasma under otherwise similar condi- 
tions. The results are illustrated using the case of the 
30<i state. When only 30d atoms are created, almost no 
electrons are detected in Gil or GI2. However, when 
the plasma is created, some electrons appear in GI2 but 
the electron number in Gil is unchanged. This indicates 
that more free electrons are formed when the Rydberg 
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laser is present. These electrons must come from the Ry- 
dberg atoms, indicating that the plasma accelerates the 
ionization of the Rydberg atoms. 
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FIG. 1: (Color online) Electron signal GI2 monitored at ti — 
4.7 /is after plasma creation when scanning the Rydberg laser. 
Upper curve: plasma present. Lower curve: plasma absent. 
In the inset: lower trace: a schematic view of the voltage 
between the two grids, the middle trace: MCP signal (the first 
peak is noise from the laser shot), the upper trace: integrator 
gates. 



Looking at Rydberg ionization efficiency (GI3 signal) 
as a function of plasma density and binding energy of 
the Rydberg state leads to the more quantitative results 
shown in Fig. [3 The denser the plasma, the more effi- 
cient the ionization. Furthermore, Rydberg atoms which 
do not ionize spontaneously can be ionized by the plasma. 
Long plasma-Rydberg interaction times favor the Ryd- 
berg ionization process. For instance, we were able to 
ionize Rydberg states with the plasma down to the dye 
limit of n = 19 at delays of 10 — 20 /is. The maximum 
effective evolution time is 20 /is, since after this delay 
plasma expansion leads to too small a density to have ef- 
ficient ionization. It is also worthwhile noting that simi- 
lar results are obtained when the plasma is much smaller 
in size (~ 0.1 mm diameter) than the Rydberg sample 
(~ 0.6 mm diameter). The Rydberg sample is still com- 
pletely ionized (for long interaction time), presumably as 
a consequence of the plasma expansion. 

All of these observations agree with the model based 
on avalanche ionization due to collisions between plasma 
electrons and Rydberg atoms. Some further observations 
help to understand the role of the other processes men- 
tioned in the introduction. First, the effect of continuum 
lowering, which reflects the fact the zero of energy of the 
isolated atoms is shifted by long-range Coulomb interac- 
tion with neighbors in the presence of the plasma, will 
ionize only n > 40 at short delays for the highest den- 
sity we could achieve ^j- This is in contradiction with 
the ionization of the 19d state we observed after 10 /is. 
Secondly, the efficiency of plasma-induced Rydberg ion- 
ization is strongly dependent on the plasma density but 



insensitive to Rydberg density, ruling out any possible ef- 
fect of autoionization of Rydberg pairs on the ionization 
process. Third, we have checked, up to 100 cm -1 above 
the ionization limit, that the ionization process does not 
depend strongly on the plasma electron energy. Finally, 
if we remove all the free electrons just after plasma cre- 
ation (at t\ = 0.1 /is), preliminary results on small num- 
bers of 30d Rydberg atoms show that after 10 /is, 0% are 
ionized if Pi < 15 /iJ. However, for Pi = 30 /iJ, 50% are 
ionized, whereas 100% are ionized if the electrons are not 
removed at ti — 0.1 /is. At first glance, this seems incom- 
patible with the scenario in which Rydberg ionization is 
caused by electron collisions. However, for Pi = 30 /i J 
when the free electrons are removed at t = 0.1 /is, the 
small number of electrons subsequently formed by black- 
body radiation or autoionizing atom pairs are trapped in 
the positive ion potential. We suggest that in such an 
environment, they are able to ionize a large number of 
Rydberg atoms. In addition, the creation rate for au- 
toionizing pairs may increase when electrons are present 
because they can stabilize the cloud expansion, and be- 
cause of the effect of attractive induced dipole-charge or 
dipole-dipole long range forces More quantitative 

experiments, to be published in [9j, that have been car- 
ried out at Colby College in collaboration with Labo- 
ratoire Aime Cotton and with a group at University of 
Virginia have observed redistribution to Rydberg states 
other than that initially populated by the laser. This 
seems to confirm the model in which Rydberg ionization 
is caused primarily by electron collisions because Ryd- 
berg redistribution is a necessary part of this process |7| . 

The second part of this letter is devoted to an im- 
portant question concerning energy conservation during 
the evolution of Rydberg atoms to plasma 0]. Elec- 
tron collisions with very highly excited Rydberg atoms 
gradually decrease the binding energy of the Rydberg 
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FIG. 2: (Color online) Rydberg ionization efficiency (in pur- 
cent) versus plasma laser power (Pi = 10 /iJ corresponds 
roughly to 4 x 10 5 ions) at ti = 1 /is. The GI3 remaining 
Rydberg atom signal is monitored. 
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atoms until they ionize. On the other hand, electron 
collisions with lower-lying Rydberg states increase their 
binding energy (superelastic collisions) . The interface be- 
tween these two regimes occurs roughly when the Ryd- 
berg binding energy equals few times the kinetic energy 
of the plasma electrons El • In the first scenario the 
Coulomb space charge is not affected and the electrons 
are probably cooled, so T e should decrease. However, 
in the second scenario the colliding electron is heated, 
or is ejected out of the plasma, and T e should increase. 
To experimentally discriminate between these two mech- 
anisms, we need to determine the electronic temperature 
of the plasma, T e . The basic idea to experimentally de- 
termine the temperature is to have instantaneous picture 
of the electron energy distribution using the short voltage 
pulse V\ to extract electrons. In order to adequately ac- 
celerate the electrons when the pulse amplitude is small, 
its duration must be longer than 50 ns. This is close to 
the electron thermalisation time and we probably also ex- 
tract some rethermalized electrons. A similar technique 
has been used in neutral atom Bose Einstein Condensa- 
tion experiments in a static external potential flfl ] . Our 
case, electrons in the plasma, is more complex due to 
the fact our potential <f> (sum of the electronic, ionic and 
external potential) depends on the number of particles 
trapped via Poisson's equation. Our experiment is sim- 
ilar to "runaway electron" experiments, except that we 
have an inhomogeneous plasma, which leads to theoreti- 
cal complexities [ljj- The number of electrons ejected by 
the voltage V\ is plotted in Fig. both for the plasma 
only, and for plasma plus Rydberg sample. Our lasers 
are not strongly focused so the cylindrical cloud is well 
approximated by a spherical gaussian symmetry for the 
ion density rii(r,t) — rP i e~ r ^ 2a For the data pre- 
sented in Fig. the sample is still gaussian (see Q) with 
a = a(t = 0) fa 250/xm. 

We first determine n? = ,„ ^W-s using the threshold 

value Ef 1 (see Fig. that is necessary to remove all the 
free electrons. Indeed, assuming there are always a few 
electrons with zero velocity in the ion potential well, the 
maximum electric field created by the ionic space charge 
is exactly 
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Preliminary comparison with ion numbers coming from 
an ionic detection on the MCP detector are in agreement 
with this measurement. 

We then determine the plasma state, immediately be- 
fore we apply V% (time delay t\), using the Poisson- 
(Landau)-Fokker-Planck (FP) kinetic equation for elec- 
trons The cold and heavy ions lead to a negligi- 
ble contribution in the electrons FP equation which is 
then is formally identical to the one governing globular 
cluster stars dynamics ^J|2(J) except that, in the non 
collisional part, there is a repulsive electron-electron po- 
tential versus the attractive gravitational field. The sys- 
tem reaches a quasi-equilibrium within a few times the 



electron-electron thermal momentum relaxation time r ee , 
which is typically few tens of nanoseconds. The main 
result obtained with our analogy with globular clusters 
is then that the phase-space density function / for this 
quasi-equilibrium is close to a (Michie-)King type distri- 
bution |2l| : 
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where ks is Boltzman's constant, E = q e (f> + ^m e v 2 and 
E t is the potential energy at infinity |22(. This defines 
what we call the uniform (r independent) King's elec- 
tronic temperature T e , which is not the same as the in- 
homogeneous velocity average temperature. Fitting the 
Monte Carlo simulations of Figs. 1 and 2 of Ref. con- 
firms that such King's type distribution is a better choice 
than a Maxwcll-Boltzmann distribution. 
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FIG. 3: (Color online) Number of electrons ejected (Gil) 
when varying the voltage Vi for the plasma alone, and for 
45d Rydberg atoms embedded in the plasma after t\ — 1 /is 
of evolution time. The dashed lines are guide for the eye to 
interpret EqQ Each data have been calibrated to the multi- 
ple shots average electron number N e , i.e. the signal shown 
is JV e xGIl/(GIl+GI2). 



The plasma temperature determination is then based 
on the solution of the self-consistent Poisson equation: 



i|_( rfh ( rit )) = _g_ 
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where rj t (r, t) 
n e (r,t) is: 
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and the electron density 
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where the proportionality factor is straightforwardly 
linked to the electron density at the cloud center, and 
to r](t) = rit(t, r — 0). More precisely a and n° are known 
from Eq. ^and we add the conditions: (i) n e (r) should be 
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non-zero at infinity 22], because even for small Vi values 
some electrons are removed, and (ii) the total number of 
calculated electrons N e (t) — J +O ° n e (r, t)47rr 2 dr should 
reproduce the observed number N e . This determines the 
values for the two remaining parameters w? e and T e in 
function of the unknown r\ parameter. The non analyt- 
ical final results may be approximated by the following 
intuitive results: the trapping depth ry^e^e has to be 
roughly equal to the trap depth, calculated assuming a 
gaussian shape for ions and electrons: v / 2q e /(y / 7r47r£:o) 
times (iVj — N e )/a 0. The experimental results (see Fig. 
[3J indicates that Ni — N e is not strongly affected by the 
Rydberg ionization. This might be an indication that the 
superelastic collisions, which would lead to electron evap- 
oration, are not a dominant process. This is also an indi- 
cation that the sample is not strongly heated when Ryd- 
berg atoms are ionized. We cannot give a more definitive 
answer because the plasma temperature is mainly given 
by the ratio Ni—N e divided by the unknown parameter r\. 
Determination of r\ from the shape of the curves in Fig. [3] 
is beyond the scope of this article, but this problem may 
probably be tackled using the collisionless Poisson-Vlasov 
"violent relaxation" theory. Nevertheless, 4<?7<12isa 
reasonable choice, and for plasma alone, analysis of Fig. 
2 of Ref. [7j indicates that r\ increases with time and that 
r\ w 10 for our experiment with t = 1 /is (see also Fig. 10 
of Ref. Without any simulation it seems difficult 

to know r\ when Rydberg atoms are present in plasma. 
Adding to the fact that a and N$ — N e are not strongly 
affected by Rydberg ionization. We could conclude that 
the temperature should not increase or decrease by more 
than a factor 5 when Rydberg atoms are added to plasma. 

To conclude, we have studied the ionization of excited 
atoms in an ultracold plasma. We have observed the fast 
ionization of Rydberg atoms embedded in the plasma. 
Deeply bound Rydberg atoms (n = 19) are ionized us- 



ing a dense plasma after some 10 — 20 fj,s interaction time. 
The only mechanism able to explain a Rydberg ionization 
efficiency of close to 100% is Rydberg-electron collisions 
from electrons trapped in the ion space charge poten- 
tial. Nevertheless the role of long-range forces cannot 
be ruled out and should be checked with further experi- 
ments. Finally, in analogy with globular cluster dynam- 
ics, we found that the plasma is always in a Michie-King 
quasi-equilibrium distribution, this will be further dis- 
cussed in a subsequent paper. T e is mainly given by 
Ni — N e which is experimentally determined by a forced 
fast electron extraction. These preliminary experiments, 
with low a temperature plasma, indicate that Rydberg 
atoms cannot drastically change the plasma temperature 
so this technique does not provide a path to reach the 
strongly coupled plasma regime. Measuring the electron 
evaporation rate with a static field or a voltage ramp are 
natural evolutions of the theory. Similar ideas have been 
published recently, |24( but based on a Maxwellian dis- 
tribution which we believe not to be appropriate in this 
system. We hope the analogy with globular clusters will 
be useful for the future of ultra-cold plasma physics and 
will stimulate links with the astrophysics community. 
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